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Abstract 

We present the Seyfert and star formation Activity in the Far-InfraRed (SAFIR) project, a 
smah (15. Ih) Herschel guaranteed time proposal performing PACS and SPIRE imaging of 
a smah sample of nearby Seyfert galaxies. This project is aimed at studying the physical 
nature of the nuclear IR emission by means of multi-component spectral energy distribution 
(SED) fitting and the star formation properties of AGN hosts, as traced by cold dust. We 
summarize the results achieved so far and outline the on-going work. 

1 Introduction 

The infrared (IR) range provides a unique probe into the coeval AGN and starburst phenom- 
ena. Most of the mid-IR (MIR) and far-IR (FIR) emission comes from dust, being of thermal 
origin. According to the unified model, the central engine and broad-line region (BLR) are 
obscured by a thick dust torus. The grains absorb the UV/optical photons and re-radiate in 
the IR. The dusty torus emission peaks in the MIR (7-30 /xm), being well characterized by 
existing facilities (eg. Spitzer , T-ReCS). It extends to the FIR, where the star formation 
(SF) related emission becomes dominant. So far this has been poorly constrained, due to the 
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limited spatial resolution and spectral coverage. According to [7], the MIR and FIR SED of 
Seyfert galaxies can be explained solely through thermal re-radiation of higher energy pho- 
tons by dust; this thermal emission is made up of three different components: (a) warm dust 
(120-170 K), heated by the AGN; (b) cold dust (40-70 K), heated by star formation, and (c) 
very cold dust (15-25 K), typical of dust heated by the general interstellar radiation field. 
Herschel [5j offers a unique window to study the FIR/sub-mm emission from nearby galaxies: 
on the one hand, the PACS photometer permits imaging at 70, 100 and 160 fim bands 
with unprecedented spatial resolution (^5.5arcsec at 70 /xm), thus allowing to characterize 
the SED minimizing the contamination from the host galaxy; and on the other, the SPIRE 
photometer allows to sample the formerly unexplored region in three bands at 250, 350 and 
500 fim with a relatively high spatial resolution, probing the cold and very cold dust compo- 
nents across the galaxy and even in the nuclear and circum-nuclear region. These data can 
be used to build SEDs samphng the peak and Ray leigh- Jeans tail of the thermal cold/very 
cold dust emission, allowing to derive masses and temperatures. The star formation rate 
(SFR) can be also derived. The spatial sampling also allows distinguishing different regions 
within the galaxies (e.g. nucleus, arm, inter-arm regions). The two main subjects of study 
within SAFIR are: (a) The dusty torus: the current models fall in two main categories: 
smooth and clumpy distributions. When applying any AGN torus model to the actual SEDs 
of Seyfert galaxies, three ingredients should be considered: AGN, starburst and host galaxy. 
The starburst contribution is mainly determined by the FIR data points. Moreover, the torus 
and starburst emissions overlap smoothly in the FIR. Therefore, the availability of FIR data 
with high spatial resolution and wide spectral coverage is of crucial importance to break the 
model degeneracy, (b) Study of star formation and nuclear activity: the interrelationship 
between star formation and accretion onto massive black holes is crucial to understanding 
galaxy formation and evolution. Herschel provides a unique means to map the star formation 
activity across the host galaxies of AGN, allowing to to determine the location and extent of 
the (obscured) star forming regions and the variation of the SFR with the radial distance, 
and permitting an accurate mapping of the the dust composition and temperature. 



2 Sample of galaxies and technical implementation 

The SAFIR sample comprises 18 nearby galaxies that were selected in such a way that 
different nuclear classes (Seyfert 1.x & Seyfert 2) were represented. Ten objects are barred 
spirals/lenticular galaxies. Five objects are peculiar/interacting systems. Four galaxies are 
confirmed Luminous or Ultra-luminous IR galaxies (LIRG/ULIRG). Only objects with either 
high-resolution ground-based MIR or Spitzer data were chosen. In addition, all the targets 
have optical, NIR, X-ray and radio data. Therefore, a complete multi-wavelength SED can 
be constructed for all targets. 

The observations were carried out using PACS and SPIRE scan map modes. The 
mapped areas for all targets were adjusted in such a way that the host galaxy and a back- 
ground region fit within the surveyed area. The la sensitivities achieved by the PACS pho- 
tometer were approximately 3.6, 3.9 and 3.9 mJy/beam at 70, 100 and 160 /xm, respectively. 
For the SPIRE photometer, la sensitivities just above the extragalactic confusion limit were 



Sanchez-Portal et al. 



3 



achieved: 5.5, 7.6 and 6.4 mJy/beam at 250, 350 and 500 /xm. These sensitivities allowed 
mapping, not only the nuclear and circum-nuclear regions but also large areas within the 
spiral discs. 

3 Results 

In this section we present some already published results from the SAFIR project O [21 [1] 
and on- going work. 

3.1 NGC 3081 

NGC 3081 is an early-type barred spiral galaxy ((R)SABO/a(r)), harboring a Seyfert 2 nu- 
cleus. A series of star-forming nested ring-like features are observed: nuclear (rl, 2.3 kpc), 
inner (r2, 11 kpc), outer (26.9 kpc). Our study combines Herschel FIR data with high- 
resolution ground-based NIR/MIR data [9J. The r2 ring is clearly resolved in our images up 
to 250 /xm. The nuclear SED was computed with unresolved FIR fluxes (r < 1.7 kpc) plus NIR 
and MIR data as shown in Fig. [l] (left panel). The torus was modeled using a clumpy model 
[1]. We wanted to test how the inclusion of FIR data affects the torus parameters, specially 
the torus size. We found that the inclusion of FIR data results in a notable increase of the 
torus outer radius: R^ = 4^^ pc vs. R^ = 0.7 =b 0.3 pc obtained using only NIR and MIR data 
[To] . The radial distribution of clouds (as given by the clumpy model q parameter) also flat- 
tens when FIR data are included in the modeling (q = 0.2 vs q=2.3). The remaining fitting 
parameters (width of the angular distribution, number of clouds, inclination angle and cloud 
optical depth) are in agreement with those obtained without the inclusion of the FIR range. 
At larger scales (1.7kpc < r < 5.4kpc), the FIR emission can be reproduced by cold dust at 
T = 28±1K (assuming a modified blackbody with emissivity /3 = 2), heated by young stars 
within the galaxy disc (likely located at rl). Finally, the FIR emission from the outer part of 
the galaxy is compatible with dust heated by the interstellar radiation field (T = 19=b3K ). 

3.2 Mrk 938 

The morphologically peculiar galaxy Mrk 938 has been proposed to be the remnant of a 
gas-rich merger of two unequal mass galaxies [TT]. According to its IR luminosity is a LIRG. 
The galaxy hosts a Seyfert 2 AGN and is known for a vigorous starburst activity. We have 
performed a multiwavelength study, including X-ray, NIR, MIR and Herschel FIR data [2]. 
The decomposition of the MIR Spitzer /IRS spectrum shows that the AGN bolometric con- 
tribution to the MIR and total IR luminosity is small [L5o/(AGN)/L/i^^0.02], which agrees 
with previous estimations. We have characterized the physical nature of its strong IR emis- 
sion. The MIPS 24 /xm and PACS 70 /xm images reveal that the bulk of the star formation 
activity is located in a compact, obscured region of < 2 kpc. We have used Herschel imaging 
data for the first time to constrain the cold dust emission with unprecedented accuracy. We 
have fitted the integrated IR SED (see Fig. [1} right panel) and derived the properties of the 
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Figure 1: Left: IR SED of NGC 3081 from 9 . Solid and dashed lines are the best fit to the data (MAP) 
and the model described by the median of the posteriors, respectively. The shaded region indicates the range 
of models compatible with the observations at the 2a level. The SPIRE nuclear fluxes include emission from 
the rl ring, and therefore have been considered as upper limits in the fit. Right: IR SED of Mrk 938 from [2], 
including modeling to the MIR and FIR data using two modified blackbodies (solid black line) , with T^^ = 67 K 
and Tc = 35K (dashed lines). The grey dot-dashed line is the fitted type 2 AGN template. 

dust. We have found that the MIR to FIR spectrum can be properly modeled with two mod- 
ified blackbodies, with fixed emissivity /3 = 2 and best-fitting temperatures of T-^^ = 67K and 
Tc = 35 K for the warm and cold dust dust components, respectively. The FIR-only spectrum 
can be fitted to a single modified blackbody with ^5 = 2 and T = 36.5 K. Using this value and 
the SPIRE 250 fim flux, we have derived the dust mass using the expression from [6J assuming 
an absorption coefficient /^250/im = 4.99cm^g~-'^. The obtained value, M^ust^^x 10'' M©, is 
consistent with values derived for local ULIRGs and other IR-bright galaxies. 

3.3 NGC 1365 

NGC 1365 (Fig. [2} top right) is a supergiant barred spiral galaxy (SB(s)b). It is a nearby 
(18.6 Mpc) LIRG harboring a Seyfert 1.5 type nucleus. The inner Linblad resonance (ILR) 
region of the galaxy contains a powerful nuclear starburst ring with an approximate diameter 
of 2kpc. We have probed the nuclear and circumnuclear activity of this galaxy in the IR 
pLj. The strong star formation activity in the ring is resolved by the Herschel /PACS imaging 
data that shows some substructures (super star clusters), as well as by the Spitzer 24 /im 
continuum emission, [Ne ll]12.81 /xm line emission, and 6.2 and 11.3/xm PAH emission. The 
active galactic nucleus (AGN) is the brightest source in the central regions up to A ^ 24 /xm, 
but it becomes increasingly fainter in the FIR when compared to the emission originating in 
the IR clusters located in the ring. We modeled the AGN unresolved IR emission with the 
clumpy torus models and estimated that the AGN contributes only to a small fraction (^ 5%) 
of the IR emission produced in the inner ^ 5 kpc. The torus size, estimated using the clumpy 
torus model [4J is ^ 5 pc. We fitted the non-AGN 24-500 /im spectral energy distribution of 
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the ILR region and found that the dust temperatures and mass are similar to those of other 
nuclear and circumnuclear starburst regions. Finally, we show that the comparison of the IR 
SFR with that derived from Ha indicates that ^ 85% of the on-going star formation within 
the ILR is taking place in dust-obscured regions. 

4 On-going work 

We have started a study of the dust properties of spatially well-resolved AGN hosts (Sanchez- 
Portal, Castillo- Fraile et al. in preparation). Within our sample, four galaxies (NGC 1365, 
NGC 4258, NGC 1566 and NGC 5728) have an apparent size large enough to allow a detailed 
analysis of the spatial dust properties, notably its temperature and mass that can be di- 
rectly compared with the star formation characteristics. For these objects, we are exploiting 
the spatial resolution of our observations to produce dust temperature maps. Assuming an 
optically thin emission, the flux density can be expressed as fi^ oc B{iy^Tdust) where /3 is 
the dust emissivity. As already stated, several dust components with different temperature 
should be generally considered, so the flux density can be expressed as: 



where N^ are the normalization constants and are the temperatures of the different com- 
ponents. The procedure devised to generate temperature maps includes the following steps: 
after a standard reduction procedure (see |lj for a description), the PACS 70, 100 and 160 /xm 
and SPIRE 250 and 350 /xm maps have been convolved to the resolution of the SPIRE 
500 fim images and resampled to the largest pixel size (that of the SPIRE 500 fim maps, 
set to 14arcsec). The images have been spatially registered and used as input to an IDL 
procedure that performs a least-squares fit to either one or two dust components (n= 1, 2) at 
each pixel. In the maps shown in Fig. [2] we have used a single temperature component with 
a fixed emissivity /3 = 2 to create the temperature maps of NGC 1365 (top) and NGC 1566 
(bottom). The latter is a bright (LC II-III), nearby (11.83 Mpc) SAB(s)bc spiral galaxy har- 
boring a Seyfert type 1.5 nucleus. By comparison with the morphology of 70 /xm and optical 
Ha images, it can be observed how closely the temperature maps follow the star formation 
regions, with the highest temperatures corresponding to areas of high SF activity. 

5 Conclusions 

The high spatial resolution Herschel PACS & SPIRE observations are demonstrating the 
importance of the FIR to improve the quality and reliability of the AGN torus fits. Parameters 
as important as the torus radius and cloud radial distribution have a strong dependency of this 
spectral range. Moreover, the FIR data are crucial to characterise the starburst contribution 
and to constrain the dust properties. The quality of the Herschel data is allowing us to 
study the spatial distribution of dust within the galaxies, thus permitting to characterize the 
variation of dust properties (e.g. temperature) and SFR with the nuclear distance. 




n 



(1) 



6 



SAFIR: testing AGN and star formation activity. 



D 



llormaliJafion 



Figure 2: Temperature maps of NGC 1365 (top) and NGC 1566 (bottom). Average temperatures range 
from ^ 17-18 K in the inter-arm regions to T ^23-24 K in the bright spots within the spiral arms. The highest 
average dust temperatures are observed in the central region of NGC 1365 with T^26K. 
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